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We have studied the adsorption of NO on small Rh clusters, containing one to five atoms, using 
density functional theory in both spin-polarized and non-spin-polarized forms. We find that NO 
bonds more strongly to Rh clusters than it does to Rh(100) or Rh(fll); however, it also quenches 
the magnetism of the clusters. This (local) effect results in reducing the magnitude of the adsorption 
energy, and also washes out the clear size-dependent trend observed in the non-magnetic case. Our 
results illustrate the competition present between the tendencies to bond and to magnetize, in small 



clusters. 
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I. INTRODUCTION 

The catalytic dissociation of NO is one of the important rate-limiting steps in the conversion of automobile exhaust 
gases into less undesirable products. Many precious metals are known to facilitate the catalytic dissociation of NO; of 
these, Rh appears to be the most efficient catalyst, i.e., the energetic barrier that must be overcome before the N-0 
bond can be broken has been found to be the lowest on Rh surfaces.^ In general, such barriers tend to decrease on 
decreasing the coordination of the metal (catalyst) atoms. For example, experiment J 2 * 3 * 4 ! have shown that the barrier 
for NO dissociation is lowered from 0.67 eV to 0.38-0.46 eV on going from the Rh(lll) surface, where surface atoms 
have ninefold coordination, to the Rh(100) surface, where atoms at the surface have eightfold coordination. Similarly, 
Ertl and collaborators have shown that NO molecules on the Ru(0001) surface dissociate preferentially at step edges, 
where atoms are less coordinated than on flat terraces Pin confirmation of their experimental observation, Hammer 
has showrP that the dissociation barrier for NO is lowered drastically, from 1.28 eV to 0.15 eV, upon moving from 
terrace atoms to step edge atoms. (However, there are exceptions to this trend, e.g., stepped Rh(100) surfaces, despite 
possessing lower coordination than flat Rh(100), display a higher barrier for NO reduction.-) This raises the question 
of whether small Rh nanoparticles, which may be expected to have a very high surface-to- volume ratio, and thus a 
high proportion of under-coordinated sites, may lower the dissociation barrier even further, relative to the value on 
Rh(100). 

While Rh clusters thus constitute promising candidates for the nanocatalytic dissociation of NO, there is however 
another factor to be considered: small Rh clusters are unusual in that they are magnetic, even though bulk rhodium 
is non-magnetic EHOD At first sight, magnetism might be expected to further enhance the catalytic activity of small Rh 
clusters: due to the familiar phenomenon that magnetism tends to increase interatomic distances, magnetic clusters 
will possess a lower effective coordination, and thus conceivably be better catalysts, than non-magnetic ones. However, 
as we will show below, the situation is somewhat more complicated than this. Thus, while the presence of magnetism 
makes issues regarding Rh nanocatalysis more complex, these systems offer a good opportunity for studying the 
interplay between magnetism and coordination number in the efficacy of magnetic catalysts. 

As a first step towards addressing such issues, in this paper we present results from a density functional theory 
(DFT) study of the adsorption of NO on very small Rh clusters, where the number of atoms n < 5. These Rh„NO 
complexes constitute the initial state for the rate-limiting step in the catalytic reduction of NO. Moreover, the strength 



of the adsorption may be expected to give some indications about how easy it is to dissociate NO on the Rh„ cluster, 
the argument being that a greater adsorption energy would indicate stronger bonding between the metal atoms and 
NO, which in turn may signal a weakening of the N-0 bond, thus facilitating dissociation. 

The main motivation of our study is to see how size affects the bonding ability of Rh clusters. Recent work on 
clusters has shown that the study of the size-dependence of their properties falls naturally into two size regimes. At 
larger sizes (hundreds to thousands of atoms) , the properties evolve smoothly towards those of the bulk; whether they 
approach the bulk properties from above or from below, depends on the property under examination, as well as, in 
some cases, the element under consideration. However, at smaller cluster sizes, one frequently observes trends that 
are non-monotonic, so that a specific size of cluster may possess desired properties, while adding or removing even 
a single atom may alter the properties drastically. Such oscillations can, however, make more apparent the factors 
governing the behavior of clusters, thus enabling one to gain insight into even the properties of larger clusters. 

As an example of a similar system whose nanocatalytic properties have attracted recent attention, we mention 
Au clusters (with 20 or less atoms) supported on magnesia, which have been found, both experimentally and 
theoretically' 10 ' 11 ' 12 ' 13 ' 14 ' to facilitate the oxidation of CO, even though bulk Au and Au surfaces are, famously, chem- 
ically inert. Similarly, the reaction rates of the hydrogenation of toluene^ and the cyclomerization of acetylene^ on 
various clusters (in some cases supported on oxide substrates) have also been found to vary sensitively as a function 
of cluster size. 

The rest of this paper is structured as follows: In Section [TX| we summarize previous experimental and theoretical 
work, both on bare Rh clusters and Rh„NO complexes. In Section ITTT] we present some of the technical details of our 



calculational method. In Section IV we present our results, first on the bare Rh clusters, and then on the clusters 
with NO adsorbed on them. In both cases, we present separately the results of spin-polarized and non-spin-polarized 
calculations. Finally, in Section [V] we summarize and analyze our results and discuss some of their implications. 

II. PREVIOUS WORK ON RH„ AND NO-RH„ 

The possibility that Rh clusters (unlike bulk Rh) may be magnetic was first suggested by the DFT calculations 
of Reddy et aZP In their work, thirteen-atom clusters of Pd, Rh and Ru were all found to have non-zero magnetic 
moments, with the largest value (of 21 /j,b, corresponding to 1.61 \xb per atom) being displayed by RI113. They 
attributed this to the reduced dimensionality as well as the high symmetry of the icosahedral thirteen-atom clusters. 
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Subsequently, other authors have found other geometries for RI113 that were lower in energy, with a lesser (but non- 
zero) magnetic moment! 17 * 18 * 19 * ^ Extending this work to other sizes, several groups have performed calculations to 
determine the structure, binding energies and magnetic properties of small Rh clusters using DFT as well as other 
quantum chemical techniques! 17 * 21 * 22 * 23 * 24 * 25 * 26 * 27 * There is considerable variation in the structures, binding energies 
and spin multiplicities obtained by various authors. For example, while many author j- 21 * 22 * 23 * 2 ^ predict that the 
lowest-energy structure of RI14 is non-magnetic with a tetrahedral geometry, other calculations 25 suggest that a spin 
septet state, also in a tetrahedral geometry, is favored. Similar issues exist also for RI15 and Rhg. 

Experimental studies on bare Rh clusters have been carried out by Cox et aZPSl w ho found that the clusters do 
indeed possess rather large magnetic moments, ranging from 0.3 - 1.1 (13 per atom; the magnetic moment per atom 
decreases with the size of the cluster, becoming zero in the neighborhood of n = 60. However, the smallest cluster 
studied by them corresponds to n = 8. The only experimental work on smaller bare Rh clusters that we are aware of 
is that of Gingerich et al. who studied Rh2,^and determined its binding energy and bond length. 

The adsorption of NO on a Rh dimer has been examined theoretically in a DFT study by Endou et al.^they found 
that the N-0 bond length is elongated (with respect to the gas phase) due to a back-donation mechanism whereby 
electrons are transferred from the dimer to the adsorbate. The same group also performed a comparative study of 
NO adsorption on four-atom clusters of Rh, Pd, Ag, Ir, Pt and Au, finding that the adsorption is most favored on 
Ir 4 P 

There have also been a few experimental studies of NO adsorption and dissociation on small Rh clusters. In their 
Fourier transform ion cyclotron resonance studies of reactions of NO with Rh^", Ford et alW^ observed biexponential 
kinetics, which they have interpreted in terms of structural isomerism. In another experimental study of NO decom- 
position on small, charged Rh clusters, Anderson et al: 31 have reported that, for both cationic and anionic clusters, 
the reaction rate increases smoothly with cluster size, though the reaction proceeds significantly faster on cationic 
clusters than on anionic ones. Motivated by these experimental studies, Harding et oZ.^have performed DFT studies 
of NO adsorption and dissociation on structural and spin isomers of Rh^ . Their studies suggest that the biexponential 
kinetics observed by Ford et al. is indeed due to the presence of structural isomers of Rh^J" , rather than due to different 
spin states. They found that the energy barrier for NO dissociation on Rh^, with a geometry of a trigonal prism, is 
markedly lower (0.23-0.36 eV depending on the spin state) than that reported for the same reaction on the Rh(100) 
surface (~ 0.5 eV^. 
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III. DETAILS OF AB INITIO CALCULATIONS 

All our DFT calculations have been performed using the PWscf code, which forms a part of the Quantum- 
ESPRESSO distribution.^ The Kohn-Sham equations^ were expanded in a plane-wave basis set with a cut-off 
of 30 Ry, while a larger cut-off of 216 Ry was used for the augmentation charges introduced by the ultrasoft 
pseudopotentialiP^ that we used to describe the electron-ion interactions. Since we wanted to investigate the ef- 
fects of magnetization, all our calculations were performed using both the spin-polarized (SP) and non-spin-polarized 
(NSP) versions of the Perdew-Burke-Ernzerhof (PBE^pOform of the generalized gradient approximation (GGA). Since 
the code makes use of periodic boundary conditions, the clusters were placed in a box of side 12 A; this size is large 
enough to ensure that the interaction between periodic images is negligible. Accordingly, Brillouin zone integrations 
were performed using only the T point. In order to hasten convergence to self consistency, we have used a very 
small Gaussian smearing with a width of 0.002 Ry - note that it is important that this smearing be small, since we 
have found that larger values of the smearing width can lead to errors in the magnetic moment of the lowest-energy 
configuration. 

As a reference, we have performed calculations on bulk Rh and NO in the gas phase. For the former, we obtain a 
lattice constant of 3.86 A and a bulk modulus of 254 GPa, which are in good agreement with the experimental values 
of 3.80 A^3 and 269 GPaJ^ while for the latter we obtain an N-0 bond length of 1.17 A. This too closely matches 
the experimentally determined value of 1.15 AP^ 

In agreement with previous work, we have found that the energy landscapes of the Rh„ and Rh„NO complexes 
possess a great many nearly degenerate local minima, both in coordinate space and in spin space. For this reason, 
we have made use of a very large number of starting configurations (both geometry and spin), and also performed 
some calculations where the magnetic moment was constrained. We are therefore reasonably confident that we have 
found the global-minimum structures. Structural optimization was performed using Hellmann-Feynman forces^iM2l 
and a BFGS-based algorithm^ for the minimization of energy. Moreover, no symmetry constraints were imposed 
when performing structural optimization, so as to ensure that distortions were permitted. 

For Rhi and RI12, the question of choosing initial geometric configurations is trivial. For the bare RI13 cluster, we 
tried different triangle-based geometries [equilateral (eq), isosceles (isos) and scalene] as starting configurations. For 
the initial geometries for Rh 4 and Rh 5 , we used the lowest-lying structural isomers reported in the literature: 21 22 23 



6 

the square (sq) and the tetrahedron (tet) geometry for RI14 and the triangular bipyramid (tbp) and the square pyramid 
(sqp) for RI15. Further, we made use of a variety of starting spin states; in a few cases we found that it was necessary 
to constrain the value of the magnetization in order to find some low-lying states. 

The situation becomes considerably more complex when considering NO adsorption on the Rh„ clusters, since there 
are a large number of inequivalent adsorption sites, as well as degrees of freedom corresponding to the orientation 
of the NO molecule relative to the cluster. We have therefore considered a great many starting geometries, and in 
the majority of cases we found that these relaxed to different local minima of the Rh„NO complex. The number 
of inequivalent possibilities that have to be considered increases rapidly as the size of the cluster grows, and it is 
therefore a rather challenging task to find the global minimum structure. 

IV. RESULTS 

Though it is well-established that small Rh clusters are magnetic, we have performed our calculations both per- 
mitting spin polarization (SP) and suppressing it (NSP). A comparison of the SP and NSP calculations should shed 
some light on the consequences of magnetism. Our results are grouped below into four sub-sections: (A) bare clusters 
+ SP, (B) bare clusters + NSP, (C) Rh„NO + SP, and (D) Rh„NO + NSP. 

A. Bare Rh clusters: Spin-Polarized 

For a single Rh atom, we obtain a ground state that agrees with that which is well-established, both theoretically 
and experimentally: a 4 F (4d 8 5s 1 ) state, with five spin-up (|) and three spin-down Q) electrons in the 4d orbital and 
one [ electron in the 5s orbital. This corresponds to a magnetic moment of 3 //,£,>. 

For the Rh dimer, we obtain a bond length of 2.25 A and a binding energy (BE) of 1.48 eV. The former is in excellent 
agreement with experiment!^ and previous calculations,^ while the latter closely matches the experimentally reported 
value of 1.46 eVp^as well as some theoretical values.^ We note that the calculated values for the BE reported in the 
literature vary over a range, depending on the level of theory used. We obtain a magnetic moment of 2 /is/atom, 
corresponding to a spin multiplicity (2S + 1, where S is the total spin of the cluster) of 5. This too is in agreement 
with earlier experiments and calculations. 

There is disagreement in the literature about whether the ground state of RI13 is an equilateral triangle with a spin 
multiplicity of 4pH] or an isosceles triangle with a spin multiplicity of 6p2l The lowest energy configuration found by 



us is in agreement with the latter. The next-lowest-lying isomer found by us is indeed an equilateral triangle, but 
its spin multiplicity is 6 and not 4. These two lowest-lying isomers are separated by only 0.008 eV/atom in BE. The 
equilateral triangle in the quartet spin state lies still higher, by an amount of 0.04 eV/atom. 

Low-lying isomers of RI14 have either a square or tetrahedral geometry. The lowest-energy configuration found by us 
corresponds to a tetrahedron with a spin multiplicity of 7; this is in agreement with one set of earlier calculations.^ 
We find two energetically degenerate isomers that lie higher than this by an amount of 0.05 eV/atom: a square 
geometry with a spin multiplicity of 5, and a non-magnetic tetrahedron. We note that some previous studies^ 21 * 23 ! 
have claimed that the latter configuration corresponds to the lowest-energy isomer. 

We find two degenerate lowest-energy configurations for RI15: both are square pyramids, one with spin multiplicity 
of 6, and the other with spin multiplicity of 8. In earlier work, one set of authors had found the former to be the 
lowest-lying isomer J^S while another had found the latter.^ We find that these lie lower than a triangular bipyramid 
with 25 + 1 = 8 by a small amount of 0.03 eV/atom. 

From these results, one can see that even at these very small cluster sizes, there are a large number of nearly 
degenerate spin and structural isomers, and it seems likely that several isomers will be simultaneously present upon 
experimentally preparing Rh clusters. 

As mentioned earlier, we are interested in seeing what effect the coordination number has on reactivity. The nominal 
coordination number of the lowest-lying isomer increases with the size of the cluster, having values of 0, 1, 2, 3 and 
3.2 as n, the number of atoms, is increased from 1 to 5. As expected, as the nominal coordination increases, the 
interatomic bond lengths increases. Since the effects on electronic structure of increased coordination number and 
longer bond lengths are expected to be linked and correlated, it is useful to define a quantity that simultaneously 
incorporates both effects. We therefore define the effective coordination number of the i th atom in a cluster byP^l 

N 'fA l ) = I ( p — \~> C 1 ) 

pRh\ K bulk) 

where the sum is calculated over all the other atoms j in the cluster, Rij is the distance between atoms i and 
j, p R \{R) is the computed spherical charge density distribution of an isolated Rh atom at a distance R from the 
nucleus, and Rbuik is the nearest-neighbor bond length in the bulk. In other words, N e ff(i) contains information 
about the ambient electronic density (due to the other atoms) that the atom i is embedded into. This is in the spirit 
of the embedded-atom method^ or effective-medium-theory,^ making the approximation that the density due to 
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the neighboring atoms can be approximated by the sum of the atomic densities. The average effective coordination 
number of a cluster is then given by (N e ff) — (1/n) N e ff(i). 

The filled black circles in Fig.[l]show how (N e ff) varies with n for the lowest-lying isomers. It can be seen that the 
variation is approximately linear, and that (N e fj) is significantly larger than the nominal coordination (filled black 
diamonds) . 

Our results for the bare Rh clusters are summarized in Table |T] where we have also compared our results to those 
of previous calculations, while in Fig. [2] we show the BE and structures of low-lying isomers. Note that both the 
spin multiplicity and the BE/atom of the lowest-lying isomer also increase monotonically with n, and that the energy 
difference between the two lowest-lying isomers is indeed minute in all cases. However, we note that the magnetic 
moment per atom decreases monotonically with increasing n, which is in agreement with the expectation that larger 
clusters, being more highly coordinated, should display a decreased tendency towards magnetization. The filled black 
circles in Fig. [3ja) show that the binding energy per atom varies more-or-less smoothly also with (N e fj). 

B. Bare Rh clusters: Non-Spin-Polarized 

When we repeat our calculations upon constraining the clusters to be non-magnetic, it becomes immediately obvious 
that the suppression of magnetism has a noticeable impact on structure. Our results are summarized in Table |TT| In 
agreement with general experience, one finds that the clusters contract upon performing NSP calculations. This is 
because in the SP case, magnetism (which is favored by Hund's rule) is in competition with the tendency to form 
interatomic bonds; upon suppressing the former, the latter tendency is increased, resulting in shorter interatomic 
bonds and hence larger values of (N e ff) (see the open circles in Fig. [IJ. While for the SP case, both isosceles and 
equilateral triangles constitute stable geometries for Rh 3 , in the NSP situation, the former relaxes to the latter. We 
note also that the near-degeneracy between the square and tetrahedral structures of RI14 appears to be lifted on 
suppressing spin polarization. 

The open circles in Fig.[3](a) show that the BE varies monotonically with (N e ff) also when magnetism is suppressed; 
this figure makes it clear that allowing magnetism stabilizes the clusters, while decreasing their effective coordination. 
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C. NO on Rh clusters: Spin-Polarized 

By considering a variety of starting spin states and geometries, we have found several stable configurations when 
NO is adsorbed on the Rh clusters. The structures of many of these are depicted in Figs. [4] |5]and[6j and the results 
are summarized in Table |III| where, for purposes of comparison, we have also presented our calculated values on the 
Rh(100) and Rh(lll) surfaces P 

We find that on a single Rh atom, it is most favorable for NO to adsorb in a "bent" configuration [Fig. |4ja)], 
while on the Rh dimer, the "vertical bridge" configuration [Fig. |3Jc)] is favored. On the Rh trimer, the configuration 
with lowest energy corresponds to one in which NO sits, perpendicularly, on the hollow site of an equilateral triangle 
[Fig.gk)]. Out of the many possible adsorption geometries for RI14NO depicted in Fig. [5] the one shown in Fig.|5je) is 
most optimal, with the nitrogen atom occupying the hollow site on one of the triangular faces of tetrahedral RI14, and 
the molecule being oriented perpendicular to the face. However, the vertical bridge configuration on the tetrahedron 
[Fig. [5](d)] is very nearly degenerate to this. For the RI15NO complex, the most favored geometry is that shown in 
Fig. |6jg), where NO sits slightly tilted on one of the short edges of the tpb. We note that when RI15 is instead in 
the sqp geometry, the lowest-energy configuration corresponds to Fig. [6jb). Though these two configurations differ 
by 0.14 eV in energy, it is possible that the latter configuration may be stabilized by kinetic barriers, and thus may 
be observed in experiments. 

Note that on RI12 and RI15, the most favored adsorption site is a bridge site, while on RI13 and RI14, it is a hollow site. 
We attempted to use our results to formulate guidelines for determining adsorption sites and geometries - for example, 
by examining bond lengths in the bare clusters and/or looking at the ambient electron density at the adsorption site, 
and by comparison with the Rh(100) surface (where NO adsorbs on a bridge site) and the Rh(lll) surface (where 
NO adsorbs in the hexagonal-close-packed hollow site). However, we were not able to determine any clear trends, 
and thus the problem of determining NO adsorption sites on larger clusters is likely to require a systematic trial of 
all possibilities, and thus considerable computational effort. 

We define the adsorption energy, E a( i s , by 

E a ds = EfthnNO — E Rh a — Eno, (2) 

where Ebh^no is the total energy of the Rh n NO complex, i?Rfr,o is the total energy of the lowest-lying isomer of 
the bare Rh„ cluster, and E N q is the total energy of the NO molecule in the gas phase. From the values listed in 
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Table III or the filled black squares in Figs.[3|b) and[7ja), we see that E a d s does not change significantly as a function 



of cluster size and/or (7V e //), though its magnitude is significantly higher on the clusters than on the flat Rh surfaces. 



complexes. It is clear that in all cases, the effect of NO adsorption is to lower the magnetization significantly - note 
that in the cases of RhNO and Rh 3 NO, the complex is actually found to be non-magnetic. Upon examining SP 
charge densities, we find that the magnetism is quenched most strongly for the Rh atoms that are bonded to the NO 
molecule; the magnetization in the vicinity of other atoms remains either essentially unchanged or, in a few cases, 
actually increases. As an example, in Fig. [9] we have plotted the difference between the | and [ densities, for (a) the 
lowest-energy NO-RI15 complex, and (b) the RI15 cluster, not in its equilibrium geometry but in the same structure as 
in (a). A visual inspection of the two figures shows that, in Fig. [9][a), the spin polarization is reduced significantly in 
the immediate neighborhood of the NO adsorption site. This becomes even clearer in Fig.^c), where we have plotted 
the difference between (a) and (b), i.e., the change in spin polarized charge density as a result of NO adsorption. 
In this figure, red and blue indicate an increase and decrease, respectively, in the degree of spin polarization. It is 
very clear that the five atoms of the Rh cluster fall into three groups: (i) in the two Rh atoms bonded to NO the 
spin polarization decreases (by about 34%), (ii) the next two Rh atoms show a redistribution of spin polarized charge 
density, with one set of d orbitals becoming more spin polarized while another set becomes less spin polarized; the 
overall magnetization for these two atoms does not change significantly, (iii) the remaining Rh atom, which is furthest 
away from NO, exhibits an increase in spin polarization (by about 13%); this is presumably because its bonds to the 
other Rh atoms are weakened (as they are now bonded to NO), and it is only in this one atom that the competition 
between magnetism and bonding is won by the former tendency. 



Figs. |7[b) and (c) show how the N-0 distance and the Rh-N distance vary with cluster size. The variation in 

the former is negligible, and within the limits of accuracy of our calculations. However, we note that, as is to be 
expected, the N-0 bond lengths in Rh„NO are always larger than in NO in the gas phase. The distance between the 
Rh and N atoms increases as n increases, possibly indicating a weaker bond between NO and the cluster. The slightly 
non-monotonic character of the graph in Fig. arises from the fact that the NO adsorption site varies with n: in 
some cases it is the bridge site, while in others it is instead the hollow site, resulting in slightly longer Rh-N bond 
lengths. 



In Fig. [8j we have shown how the spin multiplicity varies with n, for both the bare clusters and the Rh„NO 
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D. NO on Rh clusters: Non-Spin-Polarized 

The results presented in the three previous sections already hint that magnetism may significantly affect the 
reactivity of Rh clusters. This becomes clear when we re-do our calculations of the previous section for the NSP case, 



the results of which have been summarized in Table IV When magnetism is suppressed, we find that the magnitude 
of the adsorption energy E ac i s is significantly increased, and now varies monotonically with n, with the adsorption 
being most favored on the single Rh atom and least so on RI15 (see Fig [7ja)) . The open squares in Fig. ^b) show 
how E a ds varies with (N e ff) in the NSP case; note the monotonic dependence as well as the difference from the SP 
situation (filled black squares). It is seen that doing an NSP calculation leads to significantly higher values in the 
magnitude of the adsorption energy than those obtained upon performing SP calculations. 

Looking at the open squares in Figs. [3b) and (c), we see that, in contrast to the behavior of E a d s , the N-0 and 
Rh-N bond lengths do not appear to be very sensitive to magnetism (for RhNO and Rh 3 NO, the SP and NSP results 
are identical, since the complexes arc non-magnetic). 

V. DISCUSSION AND SUMMARY 

We list here the main findings of the previous section: (i) small Rh clusters are magnetic, as is well-established in 
the literature, (ii) the magnetic moment per atom (for the lowest-energy configuration) decreases monotonically with 
the size of the cluster, (iii) the average effective coordination number and average bond length both increase as the 
size of the cluster increases, (iv) as a result of the magnetism of the clusters, the bond lengths increase and effective 
coordination decreases, relative to what they would be if the clusters were non-magnetic (v) the adsorption energy 
E a ds does not display any clear trend as a function of n or (N e ff), in the SP case, (vi) however, E ac i s is larger (in 
magnitude) for the clusters than on Rh(100) or Rh(lll), (vii) the adsorption of NO strongly quenches (and, in some 
cases, eliminates) the magnetization of the bare clusters, (viii) this effect is local, being most prominent on the Rh 
atoms bonded to the adsorbate, (ix) repeating the calculations with spin polarization suppressed leads to still higher 
adsorption energies. 

Our finding that adsorption is weaker in the magnetic (and less effectively coordinated) case than in the non- 
magnetic one may seem, initally, to contradict the general expectation that lower coordination favors increased binding. 
However, this apparent contradiction arises from the fact that NO adsorption quenches the magnetism on the Rh 
clusters. This is similar to what has been observed for CO and NO adsorption on the magnetic Ni(llO) sur 

fac cFpTFsl 
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as well as for NO adsorption on a Rh monolayer in a hypothetical (bulk-truncated) structure.^ However, in the 
last of these cases (Ref. l49jl we note that the adsorption of NO completely removed the magnetization of the Rh 
atoms, whereas in our case this effect is only partial. For all the Rh„ and Rh„NO cases we have examined, the 
lowest-energy SP solution is lower than the lowest-energy NSP solution (except for RhNO and RI14NO, where the 
two are equal in energy, as the magnetization is totally quenched upon adsorption of NO). However, the difference 
in total energy between the NSP and SP Rh„ clusters is much larger than the difference in the total energies of the 
NSP and SP Rh„NO complexes, due to the quenching of magnetism in the latter. As a direct consequence of this, 
the SP adsorption energies are reduced with respect to the NSP adsorption energies. Our findings are consistent with 
those of Nayak et aZ.^who studied Rh„H2 clusters, and found that H2 binds much more strongly to a non-magnetic 
isomer of RI14 than to a magnetic one. 

Thus, lower coordination (which favors binding to adsorbates) also favors magnetism (which disfavors binding to 
adsorbates). This suggests that for Rh nanocatalysis, there may be an optimally effective cluster size, which is low 
enough to favor increased binding, while being high enough so that magnetism does not significantly reduce binding. 
The next step would be to see whether the greater binding to Rh clusters (relative to flat Rh surfaces) will also result 
in lower barriers for NO dissociation; work in this direction is in progress. 

We note that adsorption studies such as the ones presented here could conceivably yield three possible indicators of 
dissociation barriers: (i) the intramolecular (N-O) bond length (ii) the metal-adsorbate (Rh-N) bond length (iii) the 
adsorption energy E ac i s ■ One would expect that lowering of dissociation barriers might correlate with an increase in 
magnitude of (i) and (iii), and a decrease in (ii). Of these, it seems likely that (i) will not be a very reliable indicator, 
given that the N-0 distance is almost the same on Rh(100) and Rh(lll), even though the dissociation barriers on 
the two surfaces are different. Indeed, our results show almost no size-dependence for d,N-o, m both SP and NSP 
cases. We find that dah-N decreases as n is decreased, possibly indicating a stronger Rh-N and thus weaker N-0 
bond. However, the values of E a d s do not show such a size-dependent trend for the SP case. Thus, calculations of 
dissociation barriers are needed to resolve the issue of whether or not the barriers will vary significantly as a function 
of size. 

In conclusion, our results illustrate the fact that in small clusters there can be a competition between the tendencies 
to bond and to magnetize, with both effects being favored by reduced coordination. This interplay between bonding 
and magnetism has a significant influence on the size-dependent trends in the chemical behavior of such systems. 
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TABLE I: Comparison of our results (spin-polarized) on Rh n with previous calculation and experiments. The abbreviations for 
the different geometries of the bare clusters have been explained in Section |Hl| The numbers in parentheses against the bond 
lengths indicate the number of such bonds. For a single Rh atom in the gas phase, the spin multiplicity is 4. 



System Author Geometry (N e ff) BE Bond lengths Spin multiplicity 

(eV/atom) (A) (2S+1) 

Present - 3.61 1.48 2.25 5 

Rh 2 Ref.[23] - - 1.51 2.26 5 

Ref[2"2l - - 1.88 2.34 5 

Refill - - 1.46 2.28 5 

Present eq 4.75 2.03 2.46 6 

Present eq 4.67 2.00 2.41 4 

Present isos 4.14 2.04 2.56, 2.41 (x2) 6 

Rh 3 Ref[23] eq 1.99 2.42 4 

Refl23l isos - 1.94 2.53 (x2), 2.40 6 

Ref[22] eq - 2.35 2.45 4 

RefHU isos - 2.37 2.52 (x2), 2.48 6 

Rh 4 Present sq 5.34 2.37 2.36 5 

Present tet 5.04 2.42 2.53 7 

RefHD tet - 3.38 2.40 7 

Refl23l tet - 2.41 2.49 1 

RefHU tet - 2.91 2.50 1 

Rh 5 Present sqp 5.88 2.74 2.45 (sq base), 2.55 6 

Present sqp 5.28 2.74 2.44 (sq base), 2.62 8 

Present tbp 5.57 2.71 2.66 (x3), 2.52 (x6) 8 

RefHU sqp - 3.13 2.48 (sq base), 2.63 8 

RefH] sqp - 3.03 6 

RefUl tbp - 3.11 2.57 (x3), 2.46 (x6) 6 

Ref|2T] tbp - 2.97 8 
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TABLE II: Geometry, effective coordination number ({N e ff)), binding energy (BE) and Rh-Rh bond length (dnh-Rh) of Rh n 
(non spin-polarized). The numbers in the parentheses indicate the number of such Rh-Rh bonds present in the clusters. The 
abbreviations for the different geometries of the bare clusters have been explained in Section [nil 



System 


Geometry (N e ff) 


BE 


d.Rh-Rh 








(eV/atom) 


(A) 


Rh 2 




4.38 


0.80 


2.18 


Rh 3 


eq 


4.74 


1.71 


2.36 


Rh 4 


sq 


5.66 


2.24 


2.34 




tct 


5.93 


2.37 


2.47 


Rh 5 


tbp 


6.26 


2.54 


2.49 (x3), 2.53 (x6) 




sqp 


6.15 


2.59 


2.43 (sq base), 2.54 
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TABLE III: Geometry, NO adsorption energy (E a ds), Rh-N bond length (d,Rh-N), N-O bond length (djv-o) and spin multiplicity 
of Rh clusters with NO adsorbed on them. The lowest energy configuration for each system is given in bold text. For the sake 
of comparison we also present data on the Rh(100) and Rh(lll) surfaces. 



System 


geometry 


E a ds 


d.Rh-N 




Spin multiplicity 










eV/NO molecule 


(A) 


(A) 


(2S+1) 


RhNO 


Fig. 


4 


la) 


-3.23 


1.76 


1.20 


1 




Fig. 


4 


|b) 


-3.17 


1.75 


1.18 


1 


Rh 2 NO 


Fig. 


4 


(c) 


-3.10 


1.87 


1.21 


2 




Fig. 


4 


<i) 


-2.67 


1.83 


1.26 


2 




Fig. 


4 


» 


2.22 


-1.81 


1.18 


4 




Fig. 


4 




-2.46 


1.81 


1.18 


4 


RI13NO 


Fig. 


4 


;g) 


-2.80 


1.86 


1.22 


3 




Fig. 


4 


t) 


-2.71 


1.81 


1.19 


5 




Fig. 


4 


;o 


-2.21 


1.97 


1.31 


3 




Fig. 


4 


j) 


-2.89 


1.91, 1.99 


1.23 


3 




Fig. 


4 


» 


-2.90 


1.97 


1.23 


1 


RI14NO 


Fig. 


5 


[a) 


-2.09 


1.99 


1.21 


6 




Fig. 


5 


» 


-1.69 


2.09 


1.24 


4 




Fig. 


5 


c) 


-2.79 


1.95 


1.21 


6 




Fig. 


5 


d) 


-3.08 


1.21 


1.24 


6 




Fig. 


5 


» 


-3.09 


1.97 


1.23 


4 


RI15NO 


Fig. 


6 


!a) 


-2.57 


1.81 


1.19 


3 




Fig. 


6 




-2.67 


1.95 


1.21 


7 




Fig. 


6 


c) 


-2.50 


1.98, 2.03 


1.23 


7 




Fig. 


6 


;d) 


-2.48 


2.00, 2.02 


1.23 


7 




Fig. 


6 


» 


-2.42 


1.93 


1.21 


7 




Fig. 


6 




-2.79 


1.98, 2.03 


1.23 


7 




Fig. 


6 


is) 


-2.81 


1.94 


1.21 


7 


Rh(100) 


vertical bridge 


-2.59 


1.96 


1.20 






horizontal hollow 


-2.47 


1.98 


1.31 




Rh(lll) 


hollow hep 


-2.18 


2.10 


1.20 
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TABLE IV: Geometry, adsorption energy of NO (E a da), Rh-N bond length (dith-N) and N-O bond length (cIn-o) in Rh„NO 
in the non-spin-polarized case. The bond topology is similar to those indicated in the second column; however, note that these 
figures represent the SP and not NSP case. The bond lengths do differ in the two cases. 



System 


geometry 


E a da 














(eV/NO molecule) 


(A) 


(A) 


RhNO 


-Fig. 


4 


» 


-4.45 


1.70 


1.20 




-Fig. 


4 


lb) 


-4.39 


1.75 


1.18 


Rh 2 NO 


-Fig. 


4 


c) 


-4.64 


1.87 


1.21 


RI13NO 


-Fig. 


4 


w 


-4.22 


1.97 


1.23 


RI14NO 


-Fig. 


5 


fd) 


-3.23 


1.89 


1.22 




-Fig. 


5 


e) 


-3.07 


1.96 


1.24 


Rh 5 NO 


-Fig. 


6 


JK) 


-3.31 


1.94, 1.90 


1.22 




-Fig. 


6 


lb) 


-3.22 


1.98 


1.22 
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O NSP 

♦ nominal coordination 



°2 2.5 3 3.5 4 4.5 5 

Cluster size (n) 

FIG. 1: Effective coordination number (< N e ff >) vs. cluster size (n). The filled circles and open circles denote the < N e ff > 
of the clusters for SP and NSP cases respectively. The filled diamonds show the nominal coordination in each cluster. The 
solid and the dashed straight lines are guides to the eye for SP and NSP respectively. 
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FIG. 2: Binding energies vs. cluster size (for SP). The horizontal line corresponds to the BE per Rh atom in the bulk. The 
equilibrium structures of different clusters have been drawn with the spin multiplicity given in parentheses. 
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FIG. 3: (a) Binding energy for the bare clusters and (b) NO adsorption energy vs. effective coordination number. The filled 
and open circles show the BEs of Rh n for SP and NSP respectively. The filled squares and the open squares denote E a ds of 
Rh n NO complexes in SP and NSP cases respectively. 
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lowest 



FIG. 4: (Color online) Stable NO adsorption geometries on Rhi [(a),(b)], Rri2 [(c)-(f)] and RI13 [(g)-(k)] after geometry 
optimization (for SP cases). The Rh atoms are represented by grey spheres, N atoms by magenta spheres and oxygen by red 
spheres. The same color convention has been followed in Figs. [5] [6] and [9] The numbers in the figures are the bond lengths in 
angstroms. 
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FIG. 5: (Color online) Stable NO adsorption geometries on RI14 after geometry optimization (SP case). The numbers in the 
figures are the bond lengths in angstroms. See caption to Fig. [4] for color code. 
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cluster size (n) 



FIG. 7: (a) E ads (NO adsorption energy in eV/NO molecule), (b) N-O bond lengths (djv-o) and (c) Rh-N bond lengths 
(d/jfc-jv) of Rh n NO complexes as a function of n, the number of atoms in the cluster. The filled squares denote the SP results 
whereas the open squares denote the NSP results. 
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Cluster size (n) 

FIG. 8: Spin multiplicity of Rh„ as a function of cluster size, both before (filled circles) and after NO adsorption (filled squares). 
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FIG. 9: (Color online) Differences between | and J, charge densities (Sp a ) of (a) RI15NO and (b) RI15, but in the same geometry 
as in RI15NO. The difference in Sp s between (a) and (b) is shown in (c). The red isosurfaces in (c) denote an increase in 
magnetization and the blue isosurfaces represent a decrease in magnetization. Note the quenching of magnetism in the vicinity 
of NO in (a), which shows up more clearly in the blue lobes in (c). 



